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We consider a fractional generalization of gradient systems. We use differential 
forms and exterior derivatives of fractional orders. Examples of fractional gradient 
systems are considered. We describe the stationary states of these systems. 
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1 Introduction 

Derivatives and integrals of fractional order Pfl El EJ have found many applications in recent 
studies in physics. The interest in fractional analysis has been growing continually during 
the last few years. Fractional analysis has numerous applications: kinetic theories [H E] ; 
statistical mechanics [3 |HJ E]; dynamics in complex media [TUl HU H2 H31 EE1] : and many 
others. 

The theory of derivatives of non-integer order goes back to Leibniz, Liouville, Riemann, 
j> ! Grunwald, and Letnikov. In the last few decades many authors have pointed out that 

fractional-order models are more appropriate than integer-order models for various real ma- 
terials. Fractional derivatives provide an excellent instrument for the description of memory 
and hereditary properties of various materials and processes. This is the main advantage of 
fractional derivatives in comparison with classical integer-order models in which such effects 
are in fact neglected. The advantages of fractional derivatives become apparent in mod- 
elling mechanical and electrical properties of real materials, as well as in the description of 
rheological properties of rocks, and in many other fields. 

In this Letter, we use a fractional generalization of exterior calculus that was suggested 
in Refs. [To^ ll6j. Fractional generalizations of differential forms and exterior derivatives were 
defined in ^3] . It allows us to consider the fractional generalization of gradient dynamical 
systems [T7J [TH]. The suggested class of fractional gradient systems is a wider class than 
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the class of usual gradient dynamical systems. The gradient systems can be considered as 
special case of fractional gradient systems. 

In Section 2, a brief review of gradient systems and exterior calculus is considered to fix 
notation and provide a convenient reference. In Section 3, a definition of fractional gener- 
alization of gradient systems is suggested. In Section 4, we consider a fractional gradient 
system that cannot be considered as a gradient system. In Section 5, we prove that a dynam- 
ical system that is defined by the well-known Lorenz equations |2"3*1 |2"4] can be considered as 
a fractional gradient system. Finally, a short conclusion is given in Section 6. 

2 Gradient Systems 

In this section, a brief review of gradient systems and exterior calculus [TS] is considered to 
fix notation and provide a convenient reference. 

Gradient systems arise in dynamical systems theory ^J E3 El • They are described by 
the equation dx/dt = —grad V(x), where x e R n . In Cartesian coordinates, the gradient 
is given by grad V = eidV/dxi, where x = e;£j. Here and later, we mean the sum on the 
repeated indices % and j from 1 to n. 

Definition 1. A dynamical system that is described by the equations 

^ = F l (x) (* = l,...,n) (1) 

is called a gradient system in R n if the differential 1-form 

uj = Fi(x)dxi (2) 

is an exact form uj = —dV, where V = V(x) is a continuously differentiable function (0- 
form). 

Here, d is the exterior derivative [TH] . Let V = V(x) be a real, continuously differentiable 
function on R n . The exterior derivative of the function V is the one form dV = dxidV / dxi 
written in a coordinate chart (xi, ...,x n ). 

In mathematics ^H] , the concepts of closed form and exact form are defined for differential 
forms, by the equation duj = for a given form uj to be a closed form, and u = dh for an exact 
form. It is known, that to be exact is a sufficient condition to be closed. In abstract terms 
the question of whether this is also a necessary condition is a way of detecting topological 
information, by differential conditions. 



Let us consider the 1-forms (J2J. The formula for the exterior derivative d of (J2J) is 

1 /<9F; &F,\ , 
du = 2 fe ~ &^J ^ A **' 

where A is the wedge product. Therefore the condition for u to be closed is 

dF t _ dFj 
dxj dxi 

In this case, if V(x) is a potential function then dV = dxidV/dxi. The implication from 
'exact' to 'closed' is then a consequence of the symmetry of the second derivatives commute, 

d 2 V d 2 V 

(4) 



dxidxj dxjdxi 

If the function V = V(x) is smooth function, then the second derivative commute, and Eq. 
flU) holds. 

Proposition 1. If a smooth vector field F = ejFj(x) of system (Qp satisfies the relations 

(|3J) on a contractible open subset X of R n , then the dynamical system (OP is the gradient 

system such that 

dxi dV(x) 

dt dxi 

If the exact differential 1-form u is equal to zero (dV = 0), then we get the equation 

V(x) — C = that defines the stationary states of gradient dynamical system (J3J). Here, C 

is a constant. 

3 Fractional Generalization of Gradient Systems 

A fractional generalization of exterior calculus was suggested in [T^llfij. A fractional exterior 
derivative and fractional differential forms were defined ^H]- It allows us to consider the 
fractional generalization of gradient systems. 

If the partial derivatives in the definition of the exterior derivative d = dxid/dxi are 
allowed to assume fractional order, a fractional exterior derivative can be defined [THj by the 
equation d a = (dxi) a JD°.. Here, we use the fractional derivative D" in the Riemann-Liouville 
form P that is defined by the equation 

D° f(x) = ' *L [ X mdy (6) 



where m is the first whole number greater than or equal to a. The initial point of the 
fractional derivative [JJ is set to zero. The derivative of powers A; of a; is 

T(k + 1 — a) 

where k > 1, and a > 0. Note that the derivative of a constant C need not be zero. 

Let us consider a dynamical system that is defined by the equation dx/dt = F, on subset 
X of R n . In Cartesian coordinates, we can use Eq. (£Q), where i = l,..,n, x = x^, and 
F = -Fjej. The fractional analog of Definition 1 has the form. 

Definition 2. A dynamical system {Ip is called a fractional gradient system if the frac- 
tional differential 1-form 

uj a = F l (x)(dx l ) a (8) 

is an exact fractional form uo a = —d a V, where V = V(x) is a continuously differentiable 
function, and d a = (dxi) a JD°. is a fractional exterior derivative. 

Using the definition of the fractional exterior derivative, Eq. (jSJ) can be represented as 

uj a = -d a V = -{dxi) a B^V. 

Therefore, we have Fi(x) = — D".V. 

Note that Eq. (jHj) is a fractional generalization of differential form (J2J). Obviously that 
fractional 1-form u a can be closed when the differential 1-form u> = uj\ is not closed. 

Proposition 2. // a smooth vector field F = ej-i^(ar) on a contractible open subset X of 
R n satisfies the relations 

D^-D«F,=0, (9) 

then the dynamical system (QJ) is a fractional gradient system such that 

^ = -D« V(x), (10) 

where V(x) is a continuous differentiable function and ~D^.V = —Fi . 

Proof. This proposition is a corollary of the fractional generalization of Poincare lemma 
[To] . The Poincare lemma is shown (TJ3 E3 to be true for exterior fractional derivative. 



Note that the Riemann-Liouville fractional derivative of a constant need not be zero (jJJ) , 
and we have 

B«.C = — -i -C. 

r(l-a) 

Therefore we see that constants C in the equation V(x) = C cannot define a stationary state 
for Eq. (|10J1. In order to define stationary states of fractional gradient systems, we consider 
the solutions of system of the equations D%V(x) = 0. 

Proposition 3. The stationary states of gradient system J77J) are defined by the equation 

i a-m m-l m-1 



V(x)- 



n 

i=l 






E-E^ 1 ...^(II(^) fcl ) = - ( n ) 



fci=0 fc„=0 i=l 



TTie Ck 1 ...k„ are constants and m is the first whole number greater than or equal to a. 

Proof. In order to define the stationary states of a fractional gradient system, we consider 
the solution of the equation 

D«V(aO=0. (12) 

This equation can be solved by using Eq. (jBJ). Let m be the first whole number greater than 
or equal to a, then we have the solution [30] of Eq. (J12J1 in the form 

m—l 

V(x) = \xi\ a J2 a k (x 1 ,...,x i _ 1 ,x i+1 ,...,x n )(x i ) k , (13) 

fc=0 

where the a k are functions of the other coordinates. Using Eq. (|T3*j) for i = 1, ...,n, we get 
the solution of the system of equation (|T2*|) in the form (|TT|) . 

If we consider n = 2 such that x = X\ and y = X2, we have the equations of motion for 
fractional gradient system 

^ = -D«V(x,y), ^ = -D«V(x,y). (14) 

The stationary states for Eqs. (fT%)l are defined by the equation 

to— 1 m—l 



nx.^-ixj/r 1 ^ e^v = o. 



fc=0 1=0 



The Cfcj are constants and m is the first whole number greater than or equal to a. 

The Riemann-Liouville fractional derivative has some notable disadvantages in physical 
applications such as the hyper- singular improper integral, where the order of singularity is 



higher than the dimension, and nonzero of the fractional derivative of constants, which would 
entail that dissipation does not vanish for a system in equilibrium. The desire to formulate 
initial value problems for physical systems leads to the use of Caputo fractional derivatives 
[3| rather than Riemann-Liouville fractional derivative. The Caputo fractional derivative 
EH1 El E2] is defined by 

D« f(x) ~ l r f{mKy)dy (15) 

where f^ m '(y) = d m f(y)/dy m , and m is the first whole number greater than or equal to a. 
This definition is of course more restrictive that (|5|). in that requires the absolute integrability 
of the derivative of order m. The Caputo fractional derivative first computes an ordinary 
derivative followed by a fractional integral to achieve the desire order of fractional derivative. 
The Riemann-Liouville fractional derivative is computed in the reverse order. Integration 
by part of (fT5|l will lead to 

m—l T k—a 

d: J(x) = n a j(x) - £ — --/ (fc) (o+). (16) 

t^o T ( k - a + 1) 

It is observed that the second term in Eq. (|16|) regularizes the Caputo fractional derivative 
to avoid the potentially divergence from singular integration at x — 0+. In addition, the 
Caputo fractional differentiation of a constant results in zero. 

If the Caputo fractional derivative is used instead of the Riemann-Liouville fractional 
derivative then the stationary states of fractional gradient systems are the same as those for 
the usual gradient systems (V(x) — C = 0). The Caputo formulation of fractional calculus 
can be more applicable to gradient systems than the Riemann-Liouville formulation. 

4 Examples of Fractional Gradient System 

In this section, we consider fractional gradient systems that cannot be considered as a gra- 
dient system. We prove that the class of fractional gradient systems is a wider class than 
the usual class of gradient dynamical systems. The gradient systems can be considered as 
special case of fractional gradient systems. 

Example 1. Let us consider the dynamical system that is defined by the equations 

where the right hand sides of equations (J17j) have the form 

F x = acx x ~ k + bx~\ F y = (ax + b)y~ k , (18) 
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where a^O. This system cannot be considered as a gradient dynamical system. Using 

dF x dF y _ fc 

"a "a = a?/ ^ °> 

ay ax 

we get that uj = F x dx + F y dy is not closed form 

duo = —ay~ k dx A dy. 
Note that the relation © in the form 

T) a F - T) a F = 

^y 1 X ±J X 1 V u > 

is satisfied for the system (flHjl. if a = fc and the constant c is defined by c = T(l — ct)/T(2 — a) 
Therefore, this system can be considered as a fractional gradient system with a = k and a 
linear potential function 

V(x,y) = T{l-a){ax + b). 

Example 2. Let us consider the dynamical system that is defined by Eq. (JT7j) with 

F x = an{n - l)x n ~ 2 + ck{k - l)x k - 2 y\ 
F y = bm{m - l)y m - 2 + cl(l - l)x k y 1 - 2 , 



where k ^ 1 and I ^ 1. It is easy to derive that 

ckl x k - 2 y l ~ 2 {(k - l)y -(I- l)x) ^ 0, 



oF x dFy _ ,._,, ,_,, 



dy dx 

and the differential form u = F x dx + F y dy is not closed duo ^ 0. Therefore this system is 
not a gradient dynamical system. Using the conditions (JHJ) in the form 



d 2 F x d 2 F y 
dy 2 dx 2 



D 2 F - D 2 F = x - Z-JJL = o 



we get d a u = for a = 2. As the result, we have that this system can be considered as a 
fractional gradient system with a = 2 and the potential function 

V(x, y) = ax n + by m + cx k y l . 

In the general case, the fractional gradient system cannot be considered as gradient 
system. The gradient systems can be considered as special case of fractional gradient systems 
such that a = 1. 



5 Lorenz System as a Fractional Gradient System 

In this section, we prove that dynamical systems that are defined by the well-known Lorenz 
equations [221 121] are fractional gradient systems. 



The well-known Lorenz equations |2*3*1 VZQ are defined by 

dx „ dv „ dz „ 
= F — = F — F 

dt X1 dt yi dt Z1 

where the right hand sides F x , F y and F z have the forms 

F x = a(y -x), F y = (r - z)x -y, F z = xy-bz. 

The parameters a, r and b can be equal to the following values 

a = 10, b = 8/3, r = 470/19 ~ 24.74 . 

The dynamical system which is defined by the Lorenz equations cannot be considered as 
gradient dynamical system. It is easy to see that 

dF x dF v 

- z + a — r, 



dy dx 

dF x dF z 

dz dx 

9F y dF z 

- —2x. 



-y. 



dz dy 
Therefore u = F x dx + F y dy + F z dz is not a closed 1-form and we have 

duj = — [z + a — r)dx Ady + ydx Adz + 2xdy A dz. 

For the Lorenz equations, the conditions (JOJ) are satisfied: 

r\ 2 F — D 2 F =0 D 2 F — D 2 F =0 D 2 F — D 2 F =0 

As the result, we get that the Lorenz system can be considered as a fractional gradient 
dynamical system with potential function 

V(x, y, z) = -ax 3 - -ayx 2 + - (z - r)xy 2 + -y 3 - -xyz 2 + -z 3 . (19) 

The potential (|T9|) uniquely defines the Lorenz system. Using Eq. (JTTJ), we can g e ^ that the 
stationary states of the Lorenz system are defined by a = m = 2, and the equation 

V(x, y, z) + Cqo + C x x + C y y + C z z + C xy xy + C xz xz + C yz yz = 0, (20) 



where Coo, Cx, C y , C z C xy , C xz , and C yz are the constants. The plot of these stationary 
states of Lorenz system can be derived by using computer. 

Note that the Rossler system J2S], that is defined by the equations 

can be considered as a fractional gradient system with potential function 

t r / \ / \ 2 2 S 2 / \ 3 

V[x, y,z) — -(V + z)x xy y z (x — c)z . 

\ ,y 5 ) 2 \y i 2 y 3Q y 1Q g v > 

This potential uniquely defines the Rossler system. The stationary states of the Rossler 
system are defined by Eq. (J20J). 

6 Conclusion 

Using the fractional derivatives and fractional differential forms, we consider the fractional 
generalization of gradient systems. In the general case, the fractional gradient system cannot 
be considered as gradient systems. The class of fractional gradient systems is a wider class 
than the usual class of gradient dynamical systems. The gradient systems can be considered 
as special case of fractional gradient systems. Therefore, it is possible to generalize the 
application of catastrophe and bifurcation theory from gradient to a wider class of fractional 
gradient dynamical systems. Note that the order of fractional derivative can be considered 
as an additional parameter that can leads to bifurcation. For example, fractional gradient 
system with the Ginzburg-Landau potential 

-T T / \ 4 ^2 7 

V (x) = —x H — x + ox 

has the stationary states that are defined by the equation x 1 ^ a (x 3 + a'x + V) = 0, where 
bifurcation is defined by new values of parameters 

T(5 -a) u . T(5 - a) 



a — -, r, b' = b- 

6r(3-a) 6r(2 



a) 



Note also that the some of fractional gradient systems can be non-local in coordinates, 
due to the integral in the definition of fractional derivatives. These non-local properties will 
be considered in the next paper. The fractional generalization of differential forms [T5J E3 
leads us to the following open questions: Is there a fractional analog of the homology and 
cohomology theories? Is there a connection with the non-local character of the fractional 
derivative and the topological properties of the fractional differential forms? These interest- 
ing open questions require the additional research. 
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